Equation of motion approach for the Hubbard model: 
improved decoupling scheme, charge fluctuations, and the metal-insulator transition 
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A new decoupling scheme is developed for the Hubbard model which provides a unified description 
of the spin-symmetric (paramagnetic metallic and insulating) phases as well as the broken-symmetry 
AFI phase. Independent of magnetic ordering, the scheme yields, in the lowest order, the correct 
strong-coupling bandwidth of order J (for the NN hopping model) and band gap of order U, a non- 
zero critical interaction strength (above which the band gap opens) only if same-sublattice hopping 
(e.g., NNN hopping) is also present, and the correct integrated spectral weights (1/2 per spin) 
in each band for the half-filled model. The effects of charge and spin fluctuations, including spin 
twisting due to finite spin correlation length, are investigated within a static, random approximation. 
A self-consistent evaluation of the disorder-averaged self energy within the CPA is carried out 
numerically. Fluctuations activate the hopping term at first order resulting in band-broadening, 
and the consequent decrease in the band gap with fluctuation strength is obtained for several U 
values in two and three dimensions. We find that the band gap shrinks to zero continuously, 
and subsequently the density of states iV(0) between the bands grows continuously, leading to a 
continuous metal-insulator transition. For finite doping there is transfer of spectral weight between 
the Hubbard bands, and a qualitative change in the nature of the quasiparticle band dispersion, 
with an effective doping-induced hopping strength and bandwidth of order xt. 
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I. INTRODUCTION 

Electron correlation effects in narrow d- or f-bands 
|jlj are manifested in a variety of remarkable macro- 
scopic phenomena, such as metal-insulator transition in 
transition-metal compounds, magnetism and high- 
temperature superconductivity in cuprates, Q colossal 
magnetoresistance in manganates, and non Fermi-liquid 
behavior in transition- metal oxides, j^j high-T c supercon- 
ductors, H and heavy-fermion metals. (7j Incorporating 
electronic correlation in essence, the prototypical Hub- 
bard model, theoretically studied first in some detail by 
Hubbard, p-[Tof has been therefore intensively studied in 
recent years in these contexts. While much attention has 
also been devoted to the doped region away from half- 
filling, j| motivated by the discovery of high-T c super- 
conductivity in the doped cuprates, the discussion in 
the following is mainly limited to the half-filled case, fo- 
cusing only on the magnetic and metal-insulator aspects 
of the phase diagram. 

Even at half-filling the Hubbard model in three dimen- 
sions has a fairly rich phase diagram in the interaction- 
temperature (U — T) space, consisting of antiferromag- 
netic insulating (AFI), paramagnetic insulating (PI), and 
paramagnetic metallic (PM) phases. Q In addition, in 
the absence of Fermi-surface nesting (for instance due to 
next-nearest-neighbor (NNN) hopping on a bipartite lat- 
tice) a gapless antiferromagnetic metallic (AFM) phase 
also appears possible. |l2| The widely different nature 
of the constituent phases (metallic and insulating, mag- 
netic and nonmagnetic) has made it exceedingly diffi- 
cult for the entire phase diagram to be studied within a 
single theoretical framework, and over the years various 



approaches have been developed to study different parts. 
The brief account given below of some of the recent devel- 
opments will also serve to highlight the different physical 
processes of importance in different sections of the phase 
diagram. 

The study of magnetism in correlated electron systems 
goes back quite some time, to the theory of spin-waves in 
the band model of ferromagnetic metals. [T|] The corre- 
lation effect is implicit in the spin-wave excitation, which 
can be regarded as a bound state of an electron of one 
spin with a hole of opposite spin. Propelled by the dis- 
covery of nearly two-dimensional antifcrromagnetism in 
La2Cu04, |l4j | the intense activity in recent years in low- 
dimensional magnetism and spin-fluctuation effects, and 
more recently in three-dimensional magnetism, has led 
to significant progress. 

Starting with the AFI phase, a major element in this 
progress has been the unified understanding of this phase 
in the whole U/t range within a many-body-theoretical 
framework, thus providing a bridge between the weak- 
coupling Slater spin-density-wave limit and the strong- 
coupling Heisenberg limit of localized spins. For a gen- 
eralized Hubbard model with multiple degenerate or- 
bitals per site, it has been shown diagrammatically that 
within a systematic inverse-degeneracy expansion, which 
preserves the spin-rotational symmetry and hence the 
Goldstone mode order by order, identical results are ob- 
tained in the strong coupling limit (U/t — > oo) for all 
quantum corrections, order by order, as from the linear 
spin- wave analysis of the equivalent quantum Heisenberg 
antiferromagnet (QHAF). |15| Various quantities char- 
acterizing the AFI phase have been evaluated in the 
whole U/t range, both in two and three dimensions, 
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and shown to interpolate properly between the weak 
and strong coupling limits. These include the spin-wave 
(magnon) velocity, JT^Jl^ ] magnon energy and density of 
states, fl8|,[l9| sublatticc magnetization m, |l6| , [l9| - pl| and 
the local transverse spin fluctuation magnitude (S]_) = 
+ S+S7). j2|@ The rapid suppression (with in- 
creasing U) in the contribution of particle-hole excita- 
tions across the charge gap, and the rapid rise in that of 
the low-energy collective (magnon) excitations indicates 
that an effective spin picture is approximately valid down 
to surprisingly low U values. [ p2[ 

For the three-dimensional antifcrromagnct, the li- 
tispendence of the Neel temperature Tn, which deter- 
mines the magnetic phase boundary between the or- 
dered AFI phase and the spin disordered PI phase, has 
also been quantitatively studied in the whole U/t range. 
A variety of techniques have been employed including 
functional integral formalism, mainly within the static 
approximation, 23-2(| quantum Monte Carlo methods, 
p7| , p8[ dynamical mean- field theory, [^9 30 which be- 
comes exact in the limit of infinite dimensions, ]3l| ] 
Onsager reaction field theory, |5^] self-consistent spin- 
fluctuation theory, |3^] etc. Except for the last two, all 
share the common feature of yielding a Tn which ap- 
proaches, in the strong coupling limit, the mean-ficld- 
theory (MFT) result for the equivalent spin-1/2 QHAF. 
This is due to the neglect of the long-wavelength, low- 
energy magnon excitations, which reduce Tn to nearly 
70 % of the MFT value within the spin-fluctuation the- 
ory, HI in excellent agreement with recent QMC results 
for the QHAF. (34) 

The [/-dependence of Tn has been shown to be very 
closely related to that of the maximum magnon energy 
lu u1 which sets the magnon energy scale. J33| Starting 
from the strong coupling limit and going down in inter- 
action strength, w m initially increases as J = At 2 /U, the 
exchange coupling. However, as the weak coupling limit 
is approached the falling charge gap (2A) compresses the 
magnon spectrum, and when 2A ~ u> m the energy scale 
uj m turns over and starts falling, eventually approaching 
2 A in the weak coupling limit. 

We first consider traversing up in temperature, start- 
ing in the ordered AFI phase in the intermediate coupling 
regime. With increasing temperature thermal excitation 
of the low-energy magnons enhances the transverse spin 
fluctuation (Sj_), thereby reducing the AF order param- 
eter, until at T = Tn ~ u> m , the sublattice magnetiza- 
tion (S£)a vanishes and the staggered susceptibility di- 
verges, marking the onset of the PI phase. With further 
increase in temperature beyond T/v, since no further re- 
duction in the sublattice magnetization is possible, a low- 
momentum cutoff must be introduced in the sum over 
magnon modes. ]33| The corresponding length scale sets 
an upper limit to the magnon wavelength, and therefore 
determines the spin correlation length £ such that AF or- 
der exists only upto length scale £. The PI phase is thus 
characterized by an energy gap, a spin correlation length, 



and local magnetization (Si) whose direction slowly fluc- 
tuates in space and time. 

The twisting between neighboring spins relative to the 
AF alignment increases as with temperature, and 
this twist activates the hopping term at first order, re- 
sulting in a broadening of the two bands, as in the spiral 
phase studied earlier in the context of the doped Hub- 
bard model. |]35|-^7| The electronic bandwidth in the low- 
temperature ordered phase, which is of the order of the 
exchange energy J = 4t 2 /U due to the second-order vir- 
tual hopping process, is enhanced due to this twist by 
order thus providing a mechanism for coupling be- 

tween the magnetic (dis)ordering and the electronic spec- 
trum. For T >> J, the correlation length £ « 0, and 
the bands are broadened to their maximum extent of or- 
der W/2, the free electron half-bandwidth. Finally when 
T « U, there is a crossover into a semiconducting phase. 

We now consider another sequence traversing down 
in interaction strength, keeping the temperature fixed. 
In this case charge fluctuations, which are strongly sup- 
pressed in the strong coupling limit, progressively become 
more important, and when U ~ T the probability of a site 
being doubly occupied or vacant becomes appreciable, 
and the local moment defined by mf = ((n^ — «ix) 2 ) be- 
comes insignificant. Double occupancy and vacany also 
activate the hopping term at first order, broadening the 
two bands. As the band gap shrinks with decreasing U 
and the bands broaden due to this charge-fluctuation in- 
duced hopping, the two bands eventually meet, resulting 
in the metal-insulator transtion. Thus the disappearance 
of the local moment and the transition into the metallic 
phase go hand-in-hand, both driven by enhanced charge 
fluctuations. 

An explicit demonstration of the existence of an in- 
sulating state with two separated bands, independent of 
any magnetic ordering, which can undergo a transition 
with decreasing interaction strength into a metallic phase 
was first provided by Hubbard. f||Tc|] However, the de- 
coupling scheme employed within the equation of mo- 
tion approach resulted in several drawbacks (discussed in 
more detail in section III), such as inappropriate band- 
width, violation of the Fermi-liquid theory, and a van- 
ishing critical U c for the metal-insulator phase transition 
for arbitrary hopping. Of course, magnetic ordering and 
correlations, and the associated low-energy magnetic ex- 
citations were totally neglected. 

Another program, which included the possibility of 
magnetic ordering, was based on a static random one- 
body approximation to the Hubbard model. [M The 
interaction term was replaced by an effective, self- 
consistent random potential with statistical correlation 
to simulate many-body correlation effects. Only Ising- 
like magnetic excitations were included and the low- 
energy magnon excitations were neglected. A continuous 
magnetic phase transition was obtained between the AFI 
and the PI phases, whereas two crossover temperatures 
were obtained around which, respectively, the DOS at 
Ey becomes appreciable and the local moments essen- 
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tially disappear. 

In recent years the dynamical mean field theory, which 
is exact in the limit of infinite dimensions, has emerged 
as a controlled approximation for studying the phase di- 
agram in the whole parameter space. In the limit of infi- 
nite dimensions the interaction self energy becomes local, 
and the quantum many-body problem reduces to the self- 
consistent solution of a single-site quantum problem of an 
impurity embedded in an effective medium. By extending 
this approach to two sites to account for the sublattice 
structure in the AF state, a continuous magnetic phase 
transition was obtained, as already mentioned above. In 
addition, when Tn was sufficiently suppressed (made pos- 
sible by including a frustration-inducing NNN hopping 
term), a first-order metal-insulator transition line, end- 
ing in a second-order critical point was also obtained. ]3C| ] 
Several associated features of the metal-insulator transi- 
tion in V2O3 have recently been discussed within the 
DMFT approach, such as volume discontinuities and re- 
sistivity changes, |39) and optical conductivity on both 
sides of the transition. |l0| The nature of the transition 
itself continues to be of interest, and a recent study of 
the changes in the density of states near the transition 
indicates that the Fermi liquid breaks down before the 
gap opens. 

While the AFI phase is well understood, and spin 
fluctuation effects can be studied about the broken- 
symmetry HF state even in the strong coupling limit, for 
the magnetically-disordered PI and PM phases a good 
analytical starting point at the same level is still not 
available. The decoupling approach presented here pro- 
vides this for studying the interplay of electronic spectral 
properties and AF correlations in these states. Further- 
more, the approach is equally applicable to the broken- 
symmetry state, allowing for a unified description inde- 
pendent of magnetic ordering. 

The paper is organized as follows. After introducing 
the equation of motion approach for the single-particle 
Green's function in section II, and briefly reviewing the 
Hubbard I approximation in Section III, improved de- 
coupling schemes are presented in section IV which take 
into account AF correlations. Especially appropriate for 
the spin-disordered PI phase, the final decoupling pro- 
cedure in subsection C offers several improvements over 
the Hubbard I approximation. Fluctuations are included 
in section V within a static random approximation and 
a self-consistent evaluationt of the disorder-averaged self 
energy within the CPA is carried out. Fluctuation effects 
on the charge gap and the density of states are exam- 
ined in section VI by studying the coupled, non-linear 
equations for the self-consistently determined complex 
self energy. Contact with known results for the broken- 
symmetry state is made in section VII, showing the ap- 
plicability of the decoupling scheme irrespective of mag- 
netic ordering. Deviation from perfect AF ordering and 
effects of short-ranged AF spin correlations on the elec- 
tronic spectrum are studied in section VIII. The effects 
of hole doping on the qualitative nature of the quasipar- 



ticle band dispersion, and the transfer of spectral weight 
between the Hubbard bands are briefly studied in section 
IX. Some conclusions are presented in section X. 

II. EQUATION OF MOTION APPROACH 

We consider the Hamiltonian 

H = Hp + Hi = - ts a\ a a l+Sy a + n H ni -l > C 1 ) 

iSa i 

where the kinetic energy term Hq involves the hop- 
ping matrix elements ts between sites i and i + 5, and 
Hi is the local interaction term. In momentum space 

Ho = Skcr £ k a L a kcr, where e k = - X)« elkr * denotes 
the lattice free-particle energy. 

We consider the equation of motion approach, in which 
the time derivative of an operator 0(t) = e lHt Oe~ lHt in 
the Heisenberg representation is obtained from the com- 
mutator with the Hamiltonian, idtO(t) = [0(t),H] = 
e lHt [0, H}e~ lHt , and apply it to the time-ordered one- 
particle Green's function 

Gij(t) = -i<*|Toi T (t)o} T (0)|*) = 

-i[^(t)(*|a iT (t)a] T (0)^)-e(-t)(*|a^(0)a| T (t)^)] . (2) 

Taking the time derivative and using dtOit) = S(t), the 
anticommutation property {cij-f, a ]|} — <% 5 & n d the com- 
mutator 

[flit, -H] = - *^2ts a l+ s^ + Uaifriii (3) 
s 

to determine dta^it), one obtains 

idtGijit) = 8(1)5,, -J^ts G i+S j(t) 
5 

-iU(9\TrHdt)oi t (t)at^0)\9), (4) 

which is an equation for G in terms of a higher-order 
correlation function. 

Decoupling the equation of motion at this stage by 
replacing the operator in the last term by its expec- 
tation value results in the Hartree-Fock approximation 

idtG^t) = - G i+Stj (t) + UimJGyit) , 

s 

(5) 

where the densities {rii a ) are determined self consistently. 
This approximation is a reasonable starting point in the 
broken-symmetry state, when the long-time average of 
the local density does not differ substantially from 
the short-time average. For the AF state, basic features 
such as the correlation effect and the associated electron 
localization in the strong interaction limit are actually al- 
ready contained at this level, as evidenced by the charge 
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gap of order U and the narrow AF bandwidth of order 
J. Within a spin-fluctuation approach, wherein trans- 
verse fluctuations about the HF state are systematically 
incorporated, the AF state of the half-filled Hubbard 
model has been extensively studied in recent years, in- 
cluding the effects of disorder, impurities, and vacancies. 
|||,||j2l) Even the strong coupling limit (U/t -> oo) 
is accessible within this formally weak-coupling expan- 
sion. This has been explicitly demonstrated by systemat- 
ically evaluating quantum corrections to various quanti- 
ties such as the magnon energy, AF order parameter, per- 
pendicular spin susceptibility, ground-state energy, spin 
correlation length. p"5[ 

However, in the spin disordered state this approxima- 
tion breaks down for strong interaction, when the expec- 
tation value of the local density (riij) — n/2 (the mean 
density) differs substantially from the short-time average 
(either nearly zero or nearly one). In this case the local 
magnetization strongly fluctuates between I and -I, so 
that the average local magnetization (mi) = (njf — riii ) 
vanishes, but the local moment /ij defined by /z? = (mf ) 
remains appreciable. This strong local moment arises 
from the correlation effect of an electron locally exclud- 
ing the opposite-spin electron due to the strong electronic 
repulsion. 

This provides a clue that the decoupling procedure can 
be improved by ensuring that the operators which are re- 
placed by their averages involve minimal fluctuations. In 
the following we describe several such decoupling proce- 
dures, all of which are quite distinct from that employed 
by Hubbard. J8| This approximation, commonly referred 
to as the Hubbard I, is first briefly reviewed here for com- 
parison. 



III. HUBBARD I 

In Hubbard I the next step was to continue 
with the equation of motion approach and exam- 
ine the time derivative of the last term in Eq. 
(4) involving the self-energy correction, ^ij(t) — 
—iU{$>\Tnn(t)ait(t)a , j^ (0)|*). Determining the time 
derivatives in dt[riii(t)ai-\(t)} from Eq. (3) and the com- 
mutator 



mi 



H } = -J2 ts ( a li a *+n ~ a Lsi a *i) > ( 6 ) 



and using the anticommutation property {ni^a^, a^} = 
Sijtiii, one obtains 

id t T t] {t) = US(t)(^\{n zl a th a^}\^) 
+ [7(*|T5 4 {n i; (t)a iT (t) + n a (t)a iT (t)}a) T (0)|*) 
= S(t)S ij U(n il ) + iUj2 t s (*|Tn i; (t)a i+5T (£)at T (0)|*) 

5 

- iU 2 {y\Tnl(t) ail (t)al^O)\y) + lUj^ts (*|T 



H(*k+<u(*) - 4-^(*)M*)}m*kV°)I*> • ( 7 ) 

In order to obtain a closed equation for T, the 
equations were decoupled at this stage by notic- 
ing that n?i = riii, ana - replacing the operators 

Uii(t), and {al^(t)a i+ si(t) - a]_ si (t)aii(t)} in the last 
three terms by their expectation values in a spin- 
symmetric, translationally-invariant ground state, with 

(at ; (*)a <+ «5x(t) - a i -5i{ t ) a il( t )) = °) and i n il) = Kt) = 
n/2, the average density for each spin. Fourier trans- 
forming to frequency and momentum space using 



k J 

¥(*) - E / ^ i(k 



z(k.rij — ujt) 



E 

6 



k 

E« 

k 



<5 



; j(k.r i; ,-wt) 



one obtains 



Un 



[1 + e k G(M] + UT(kiu) 



(8) 



(9) 



Substituting for T(kuj) = ^[1 + e u G(\sw)]/(u - U) from 
above into Eq. (4), and solving for G(kui) one obtains 
the Hubbard I result 

uj - 17(1 - n/2) 



G(kaj) = - 
V ; {lu ~ e k )(uj - U) - 

In the half- filled case (n = 1) this can be written as 



(10) 



G(kw) 



1 _ ik. 

1 E k 



w-R 



(i) 



uj-R 

(1,2) 



(2) 



(11) 



where £ k = V^ + 4> and E l = K £ k + U)T Eu}/2 
are the two band energies. In the strong coupling limit 
these band energies reduce to £k/2 and U + 6^/2 to low- 
est order, giving a bandwidth of W/2, where W is the 
free-particle bandwidth. Thus in going from the non- 
interacting limit to the strong-coupling limit, the band- 
width decreases only from W to W/2, whereas one ex- 
pects a correlation-narrowed bandwidth of order J for 
the NN hopping model. This is a serious drawback of 
the Hubbard I approximation in the strong correlation 
limit, and can be traced to the replacement of n,-^, a 
strongly fluctuating quantity, by its average value. An- 
other drawback pertains to the absence of a critical in- 
teraction strength for the energy gap. For a general free- 
particle energy e^, the electronic spectrum is split into 
two bands for all U. Also, in general, the integrated spec- 
tral weight in each Hubbard band is not exactly one-half 
per spin, although the total for both bands is identically 
one. This amounts to an inconsistency in the insulating 
state in that the number of states in k-space lying below 
the Fermi energy does not exactly yield the number of 
particles. 
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IV. IMPROVED DECOUPLING SCHEMES 

We adopt a different approach here which removes all 
of the drawbacks cited above. We continue with Eq. (4) 
for d t Gij and consider the second derivative 



-iU mTd t [ nil (t)a n (t)}a^(0)^) 
-iUS(t) (^\{ nil a ih a]^). 



(12) 



The derivative of Gi+s,j(t) is substituted from Eq. (4), 
the derivative dt[nn(t)aii (t)} obtained as before from 
Eqs. (3) and (6), and in the last term {n^atf , aj^} = 
SijUii- Putting all this together, and multiplying by i, 
we obtain 

i 2 d^Gij{t) = idtSityij+SitWijUimi) 

5 SS> 

+ iU^t s (*|T{n i i(t)+ni +4 i(t)}a i+ 5 T (t)at 1 .(0)|*) 

-iU 2 (^>\Tnl l (t)a n (t)a} 1 (Om+tUY / ts (*|T 

s 

{«li(*)ai+5iW-«l-5iW«iiW}«iTWa]T(°)l*) ■ ( 13 ) 

We now proceed further along three different routes by 
introducing three decoupling schemes termed A, B, and 
C in the following subsections. Schemes A and B are 
most suitable for hopping terms which only connect sites 
of opposite and same sublattices respectively. The best 
features of A and B are then combined in the final scheme 
C which is applicable for any general hopping. However, 
as we shall see, we are constrained to the consideration 
of a bipartite lattice. 



Adding -nU id t G tJ (t) + (nU/2) 2 Gij(t) to both sides 
of Eq. (13), where n = (n^ + ni{) is the total density, 
and substituting for dtGij (t) from Eq. (4), we obtain 

i 2 d 2 G lJ {t) - nU id t G l0 {t) + (nU/2) 2 G lJ {t) = id t 5(t)5 l0 
+ 5{t)[5 ij U{(n i {) - n} - ^t s S i+S ,j] + ts'Gi+s+s' ,j 

8 88' 

+ iUY / ts (*|T{n a (t) + n l+sl (t) - n}a l+5] (t) ] T (0) | *) 

8 

- iU 2 (*|T{n a (i) - n/2} 2 a lT (t)a] T (0)|*) +iU^t s 



(*\T{aUt)a i+sl (t) 



l i-8i 



(i)a a (t)H T (t)aj T (0)|*) . 



In order to decouple this equation for G involving sev- 
eral higher-order correlation functions we now replace the 
operators {riii(t) + n i+ si(t) — n}, {n»j.(t) — n/2} 2 , and 

{al^(t)ai+8 l(t) — a i-si{t) a il(t)} m the l as t three terms by 
their ground-state expectation values, again in the spin- 
symmetric, translationally-invariant ground state, with 
(riii) = ( n i}~) = n /2, the average density for each spin. 
At half-filling (n = 1), we have {n tl (t) - 1/2} 2 = (1/2) 2 , 
a c-number, so that no approximation is involved in this 
replacement. Furthermore, in the strong coupling limit, 
when double occupancy and vacancy are both prohibited, 
the total charge density operator n^ + = 1, the unit 
operator, so that the term {njj.(t) — 1/2} 2 = mf/4 is also 
related to the second moment of the local magnetization 
nii. In general, rij^(t) — 1/2 fluctuates between the two 
values nii/2 — =Fl/2, depending on whether the site is 
occupied or empty. 

The average (n^ + n^gl ~~ n ) vanishes in the param- 
agnetic state when (n^) = (m+si) = n /2- More impor- 
tantly, it also vanishes in the presence of AF correlations, 
provided the sites i and i + 5 are in opposite sublattices, 
in which case («i+5x) = ("it)- This scheme is there- 
fore most appropriate when the hopping term ts connects 
sites i and i + S of opposite sublattices, as in the NN hop- 
We take (aUt)a l+5l (t) - a\_ g ,(t)a a (t)) = 0, 



ping case. 



as in Hubbard I. This is obviously so if translational sym- 
metry is present, but also in presence of AF correlations, 
as both terms are off-diagonal in the sublattice basis, 
and of equal magnitude. After Fourier transformation 
to frequency and momentum space, and introducing a 
Hartree shift in the energy axis through the transforma- 
tion uj — nU /2 — > u), we finally obtain 



G(fcw) = 



OJ + e fe 



co 2 -(A 2 +e 2 ) 



1 



E k 



1 



+ 



E k 



w - E k uj + E k 



(15) 



where E k = ^A 2 + e{ and A 2 = U 2 ({n lL (t) - 1/2} 2 ) = 
J7 2 /4. The energy spectrum thus splits into two bands, 
independently of whether there is long-range AF order 
or not, with an energy gap of 2A = U. Furthermore, the 
bandwidth is correctly of order J in the strong coupling 
limit. 



An alternate decoupling scheme is more appropri- 
ate when the hopping term involves sites in the same 
sublattice. Adding to both sides of Eq. (13), the 
terms — nU id t Gij + (nU/2) 2 G io - nU^s^s + 
2iY,sts dtGi+sj +J2s5't$ ts ' Gi+8+8',j, and substitut- 
(14) ing for d t G(t) from Eq. (4), we obtain 



•5 



i 2 d 2 G tJ - nU id t G i:j + (nU/2) 2 G i:j 
+ 2% ^ tg d t G i+ s,j -nU^^ts Gi+gj +y]tg tg>G l+ g + g 



= id t 6(t)5ij + 5(t)[5ijU{(nn) - n} + ^^tg Si + gj\ 

s 

+ iUj2ts (*|T{n i i(t)-n <+ «i(t)}a i+aT (t)at 1 .(0)|*) 

- ill 2 (*|T{n a (t) - n/2} 2 ai i(t)a^(0)\y) + il/J^ts 

s 

WTialWcH+siit) - a|_ tfi (t)a ii (t)}a it (t)at t (0)|*) 

We notice that on the right-hand side the second-order 
hopping term (tg tgi) present in Eq. (15) has exactly 
cancelled. Again, as before, the operators are replaced 
by their expectation values. We note that we can safely 
set (rijj, — Ti»+5j.) = 0, not only in the symmetric paramag- 
netic state, but even in the presence of AF correlations, 
provided that sites i and i + 5 are in the same sublattice. 
Fourier transformation leads to 

(lu - nU/2 - e k ) 2 G{kuj) = [to - nU/2 - e k ) + A 2 G{kuj) . 

(17) 

With the same Hartree shift employed through the trans- 
formation w — nU 1 2 — > u), we obtain 



G{kuj) 



LU-e k 



1 
2 



+ 



e k + A 



scheme B is suitable when hopping is limited to the same 
sublattice. This immediately suggests that for general 
'■' hopping both schemes can be further improved upon by 
taking an intermediate measure. 

For this purpose we are constrained to the considera- 
tion of a bipartite lattice, for which the hopping terms 
can be divided into two groups, connecting sites of oppo- 
site and same sublattices. We therefore write, 



(16) 



tfj.a 



+ £k Al )aL a k < r 



(19) 



(18) 



Thus A causes the free-particle spectrum to be split into 
two bands, with no correlation effect on the bandwidth. 
The two bands overlap when the band separation 2 A = U 
is smaller than the free particle bandwidth, but a gap 
opens up beyond a certain critical interaction strength. 
Hence for a general ek, there is a non-zero critical inter- 
action strength. However, the number of states in each 
band is precisely 1/2 per spin. 



Finally, we consider a third scheme which combines the 
best features of A and B. We recall the replacement of op- 
erators by their expectation values in the decoupling pro- 
cedures, and focus on the operators {nii(t)+ni+g[(t)—n} 
and {nn(t) — n i+ si(t)} in Eqs. (14) and (16), particu- 
larly in view of AF correlations in the ground state. For 
6 connecting sites of opposite sublattices, the short-time 
average (n^ (t) + Tii+g j (t) — n) is weakly fluctuating, but 
not so when sites i and i + S arc in the same sublat- 
tice. Therefore, scheme A is suitable when hopping is 
limited to opposite sublattices. On the other hand, the 
short-time average — is weakly fluctuat- 

ing when S connects sites of the same sublattice, so that 



kcr 



where sites i + A and i + fi refer to neighbors of site i in 
the opposite and same sublattices respectively, and fur- 
thermore e kA = -£) A iAe ik - rA and e kAl = - *** eik ' r " 
are the two associated free-particle energies. 

The intermediate measure which leads to an improve- 
ment over both schemes A and B is simply effected 
by including only those additional terms of scheme B 
which correspond to the same sublattice. Thus, we add 
to both sides of Eq. (13) the terms — nU idtGij + 
{nU/2) 2 Gij - nl/£„t„ <3 i+/i j + 2 d t G l+M + 

Z)^*' V ^*+m+m'>j'' where the sum over neighbors 
(fi, u') in the last three terms are restricted to sites in 
the same sublattice. As before, we obtain after substi- 
tuting for d t G from Eq. (4), 

i 2 d 2 Gij - nU id t G tJ + {nU j2) 2 G i3 - + 2^^^ At d t G l+M 

fj- 

H mi' 

= id t 6(t)5 ij +6(t)[S ij U{{n il )-n} 

H A A, A' 

+ iUJ2^ (*|TKi(i) - n i+ ^(t)}a i+ ^(t)a^(0)\V) 
fj- 

+ tU^x (*\T{n lL (t) + n t+XL (t) - n}a 4+AT (t)a] T (0)|*) 

A 

- iU 2 (*|T{rm(t) - n/2} 2 a 4T (t)a] T (0)|*) + ill^s 

s 

(*|T{at i (t)a i+4i (t)-at_ ai (t)o ii (t)}a iT (t)at T (0)|*) . (20) 

The cancellation of the second-order hopping terms, re- 
ferred to earlier below Eq. (16), is only partial now, and 
only terms involving opposite sublattice hopping finally 
survive. 

Again replacing the operators with their expectation 
values as before, with (njj. — ni +fJ _i) = (riy + n, l+ \i —n) = 
0, we obtain after Fourier transformation 

(to - nU/2 - e^) 2 G(kuj) = u- nU/2 - e kAJ + e kA 

+ (A 2 +6 2 A )G(M (21) 
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which, after the same Hartree shift oj — nU /2 — > ui, yields 



G(kuj) 



1 

2 



H 



^A 



E- 



kA 



Ek^t 



kA 



(22) 



where i?kA = \/A 2 + e^. A is associated with the opposite 
sublattice hopping. Results of both schemes A and B are 
recovered by taking appropriate limits. For only opposite 
sublattice hopping = 0, so that Eq. (15) is obtained, 
whereas for only same sublatttice hopping £kA = 0, which 
yields Eq. (18). In general, when hopping involving same 
sublattice sites (e.g. NNN hopping) is present, a finite A 
is required for the band gap to appear, yielding a criti- 
cal interaction strength which separates the metallic and 
insulating phases. 

Quite generally, ekA, the free-particle band energy as- 
sociated with opposite-sublattice hopping, changes sign 
under the transformation k — > k + 7r, while the band 
energy ek^ T -EkA does not. Therefore, the k-space is 
split into two degenerate zones. For evaluation of the in- 
tegrated spectral weight, the k-sum therefore covers the 
whole k space. As 6kA is odd, the integrated spectral 
weight in each band is identically one-half per spin. 



V. FLUCTUATION EFFECTS 

In the previous section, the equation of motion for 
Gij(t) was decoupled by replacing the operators {mi + 
n *+Ai - n}, {n n - n i+fll } and E«{ a a (*)«»+«! (*) ~ 

a l-si(t) a il(t)} by their (vanishing) expectation val- 
ues in the spin-symmetric, translationally-invariant 
ground state. If A represents these operators, then 
{^\A(t)A(t')\^} and (^\A(t) 2 \^) are in general non- 
vanishing due to fluctuations and correlations present 
in the ground state. In the following we approximately 
account for these fluctuations by replacing the operator 
A by static, independently random terms Ai distributed 
across the lattice with (Aj) = and (Af) determined 
self-consistently from (^\A(t) 2 \ty) . 

When the above replacements are made, we see from 
Eq. (20) that the charge/potential fluctuation terms 
U (mi + n i+\i ~ n ) an d U(mi — rti+pj,) couple with the 
hopping term, leading to charge fluctuation induced hop- 
ping. On the other hand, the hopping induced fluctua- 
tion term Y,s t s( a li( t ) a i+si( t ) - a l-5i(*) a *iW) is P ur °ly 
imaginary and couples with the local interaction term 
yielding a hopping induced relaxation. 



A. Charge fluctuation induced hopping 

In Eq. (20) the operators {mi + n i+ \i — n} and 
{nil— n i+nl \ were replaced (in the decoupling procedure) 



by their expectation values, which vanish in the paramag- 
netic state at half filling (n = 1), and also when AF corre- 
lations are present. As shown below this approximation 
essentially amounts to a neglect of local charge fluctua- 
tions and is valid in the low-temperature, strong-coupling 
limit, when double occupancy and vacancy are strongly 
suppressed. In the strong coupling limit we have a spin 
picture of the ground state, and for arbitrary orientation 
of the local spin direction, the ground state expectation 
values (nil + n i+Ai) arL d ( n ii ~ remain invariant 

due to the AF correlations between neighbouring spins. 
If the spin at site i is inclined at angle 9 from the — z 
direction, then (mi + n i+ \i) = cos 2 9/2 + sin 2 9/2 = 1, 
and (nii — m+^i) = 0. Hence the operators mi + rii+Ai 
and mi — n i+ ^i essentially behave like c-numbers, and 
the above replacements are valid. Even the strong quan- 
tum spin fluctuations, which significantly reduce the AF 
order parameter in two dimensions, do not substantially 
alter this picture as the major spin-fluctuation contribu- 
tion arises from long-wavelcngth magnon modes. 

However, for weak correlation and/or high tempera- 
ture, when double occupancy is not strongly penalized, 
the local densities (mi) an d (m+5l) independently fluc- 
tuate between zero and one at short time scales, and 
therefore (mi + f^+Ai — n ) and (mi — rii+^x) strongly 
fluctuate between minus one and plus one. In the fol- 
lowing we account for these fluctuations within a static 
random approximation, together with disorder averaging 
at the CPA level. A more complete theory which includes 
the spatial and temporal correlations in the fluctuations 
will be developed later by considering the appropriate 
two-body Green's functions in Eq. (20) within system- 
atic approximations like the RPA. 

Even in the absence of charge fluctuations, (mi + 
n i+\i — n) can fluctuate if ferromagnetic correlations are 
present in the (paramagnetic) ground state. Similarly 
these fluctuations are also present if the AF correlations 
are only short-ranged, so that the neighbouring spins are 
twisted from perfect AF alignment. This is further dis- 
cussed in section VIII. 

We note that the energy scale Ut of these fluctuation 
terms is intermediate between the local U 2 scale and the 
second-order hopping scale t 2 in Eq. (14). The effect 
of these order- Ut fluctuation terms is therefore to renor- 
malize the bandwidth by order t in the strong coupling 
limit, as opposed to the order J renormalization by the 
second-order virtual hopping process. Thus the charge- 
fluctuation and spin-twisting processes activate the hop- 
ping term at first order. 

For simplicity, we consider the NN hopping case, so 
that the sites i and i + 5 are in opposite sublattices. We 
assume that U times the density fluctuation term U(mi + 
m+5l — n ) , which we denote by Vis, is a random potential 
distributed independently across the lattice. The average 
of this random potential term vanishes, but its second 
moment (V 2 S ) = U 2 ((mi + m+5l — n ) 2 ) reflects a mea- 
sure of fluctuations in the system. At half filling (n = 1), 
the long-time average yields the fluctuation magnitude 
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Sn 2 = {{riii + n i+si — n ) 2 ) = ^{ n ii n i+Si) , which has the 
following limits. In the strong coupling limit when double 
occupancy is strongly suppressed and AF correlations are 
present, we have {riiirii+si) — (iHlTHt) = 0. At the other 
extreme, in the non-interacting limit, when a site may be 
independently occupied or vacant with equal probability 
of one-half, we have 2{mini+si) = 1/2. Therefore, the 
fluctuation magnitude lies in the range < 5n 2 < 1/2. 

The random potential term Vis is treated at the CPA 
level, and the disorder-averaged Green's function is ob- 
tained in the following. After the same decoupling ap- 
proximation made earlier, we rewrite Eq. (14) in a matrix 
representation (constructed from the site basis) as 

[u 2 l - (A 2 1 + T 2 ) + £V)]G(w) = wl + T , (23) 

where the matrix T is associated with the NN hopping 
term and has matrix elements Ti t i-\-g = —t, whereas 
V has matrix elements ~Vij+s — Vi$, the random po- 
tential term. For further convenience we introduce an- 
other Green's function g(w) through the relation G(u>) = 
g(w).[wl + T], so that g(w) obeys the equation, 

[w 2 l- (A 2 l+T 2 )+*V]g(u;) = 1 . (24) 

The usual CPA treatment of the random potential 
term, involving perturbative expansion followed by dis- 
order averaging, leads to the following result at the level 
of rainbow diagrams, 

(gH> = lu 2 1 - (A 2 1 + T 2 ) - S(o;) (25) 

where, 

SH = < 2 (VgHV) (26) 

is the disorder-averaged self energy, to be determined self 
consistently. 

For the purpose of disorder averaging, we assume for 
simplicity that the random potential terms Vis are inde- 
pendently random for each pair of sites i and i + S, so 
that (Vis.Vjt') = l{SijS S s' + 5 i>j+ s'5 i+s>j ), where the sec- 
ond moment 7 = {V 2 i+S } is a measure of the potential 
disorder. This essentially amounts to a neglect of spatial 
correlations in the fluctuations. In this case there are two 
terms in the disorder self energy 

T} u {uj) = -ft 2 J2(Si+Sa+s{^)) 
s 

£ 2 M+(5 (w) = 7 < 2 (g l+ ^H) . (27) 

The second term X 2 , however, vanishes because 
(gi+s,i) = Sk'(3(k')) elk ' rs is identically zero. This fol- 
lows because (g(k')), given by the Fourier transform of 
Eq. (25), is even under the transformation k' — > k' + tv, 
while the other term e lk ■ ra is odd, and therefore the k' 
summation over the whole Brillouin zone yields the van- 
ishing result. The disorder self energy is therefore diago- 
nal in the site basis, and hence momentum independent. 



Furthermore, the sum ^ s in the equation for H 1 a{u}) 
yields z, the number of nearest neighbours, so that the 
self energy involves the combination zt 2 . These features 
of the local self energy are as in the DMFT where the hop- 
ping term is taken to be t — t*/y/z, so that zt 2 — {t*) 2 
is finite in the limit of infinite dimensions. |HJ 

Finally, after Fourier transformation of Eq. (25) and of 
{G{cj)) — (g(w)).[u;l + T], we obtain the following result 
for the full Green's function 

^ - Z^r^W) ■ (28) 

where the disorder-averaged self energy is obtained self 
consistently from 

s (") = ^E ^- (A a4)-S( M ) - (29) 

The above provides a self-consistent scheme for determin- 
ing the Green's function when fluctuations are included. 
In the following we quantitatively examine the fluctua- 
tion effects on the spectral properties, in particular the 
charge gap and the density of states. Thermal contribu- 
tion to charge and spin fluctuations brings in temperature 
into the problem, and makes the charge gap sensitive to 
temperature, leading to the possibility of a temperature- 
driven metal-insulator transition. Temporal and spatial 
correlations in fluctuations are neglected in this treat- 
ment. The present study is nonetheless useful because 
when these correlations are included by considering the 
appropriate two-body Green's functions (having the form 
E.G), the structure of the resulting equation for G will 
be similar to Eq. (28). 

B. Hopping induced relaxation 

The purely imaginary hopping induced fluctuation 

term UJ2s t ^ a \i 0H+«.l(*) - <4-5;0Hi0)) = iV l cou- 
ples diagonally to the Green's function. Therefore, Eq. 
(23) will contain, in addition, a diagonal term tiW , where 
(V')ij = &ijV(. Again assuming that V( are indepen- 
dently random, and are uncorrelated with the Vis, within 
the CPA we obtain another diagonal self energy with a 
negative coefficient which will simply renormalize the 7 
defined earlier. 



VI. CHARGE GAP AND THE 
METAL-INSULATOR TRANSITION 

A qualitative criterion for the metal-insulator transi- 
tion is obtained by considering the competition between 
the kinetic energy term (bandwidth W) which delocalizes 
the electrons making the system metallic, and the corre- 
lation term (Coulomb barrier U) which localizes the elec- 
trons making the system insulating. We will first make 
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a simple analysis of the effective charge gap when fluc- 
tuations are included to their full extent. We shall see 
that the qualitative criterion U c ~ W for the vanishing of 
the gap, and therefore for the metal-insulator transition, 
arises in a very natural way. 

We have seen earlier that when fluctuations are allowed 
to their full extent, 5n 2 = 1/2, and therefore the disorder 
strength 7 = U 2 /2. Considering the disorder self energy 
from Eq. (29) to the lowest order, we have for lo = 

^ — Er^A^ry--^ 2 , (30) 

where the result at the right is obtained by neglecting 
the free-electron energy e^/ in comparison to A = U/2, 
for the purpose of obtaining an estimate. Now from Eq. 
(28) the charge gap vanishes when |X(0)| = A 2 , which 
leads to the following estimate for the critical value U c 
for the metal-insulator transition, 

U c = V8zi? = 2V2 t* (31) 

which is of the order of the free-electron band width. 
Within the approximations used, for U > U c the energy 
gap never closes, and the system remains insulating. 

We now obtain the self-consistent solution for the self 
energy S(w) from Eq. (29). As X(w) is complex in gen- 
eral, it is convenient to introduce a complex function 

F(a, (3) = TZ(a, (3) - i/3X(a, (3) = £ 2 * 

Ij^i + ay + p Z. (e 2, +a) 2 + /32 ■ 

By equating (A 2 - lo 2 ) + E(w) = a + i(3 in Eq. (29) 
the self-consistency equation for £(w) may be recast in 
terms of a, (3\ if T = zjt 2 denotes the prefactor, Eq. (29) 
then reads — £(u;) = TF(a,f3), the real and imaginary 
components of which are 

A 2 - lo 2 = m(a, (3) + a (33) 
/3 = T(31(a,(3) . (34) 

We now discuss the solution for the above two coupled, 
non-linear equations in a and (3, where the parameter 
T = zjt 2 = zt 2 U 2 ((nn + n l+si - n) 2 ) = zt 2 U 2 Sn 2 in- 
corporates the fluctuation magnitude. We first consider 
the solution with (3 = 0, which corresponds to a van- 
ishing imaginary part of the self energy and the density 
of states, and therefore to the energy gap region. Since 
T(a,0) decreases monotonically with positive a, for a 
given T 




Q n 1 1 1 1 1 1 1 1 1 1 1 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

a/A 2 



FIG. 1. Plot of T1l{a, 0) + a vs. a, both in units of A 2 , for 
different values of the fluctuation strength Sn — 0.1 (lowest), 
0.2, 0.3, 0.4, 0.5, and 0.6 (highest). The minimum corresponds 
to a lo value from Eq. (33) which yields the band edge. As 
the minimum approaches 1 from below, the band gap shrinks 
to zero. 

there exists a physically significant solution a = a* > 
of the equation 1 = Tl(a,0). For a > a*, we have 
I(a,f3) < l(a,0) < T(a*,0), therefore Eq. (34) has 
a solution only if (3 = 0, whereas for a < a* a solu- 
tion exists only for (3^0. Thus, a* is the point where 
the self energy changes from purely real to complex, and 
therefore corresponds to the band-edge energy, as further 
discussed below. 

There is another special value a — a\ this one hav- 
ing to do with Eq. (33). The right-hand-side term 
TlZ(a, 0) + a diverges in both the limits a — > and 
a — > 00, and therefore it must have a minimum at some 
intermediate value a — a\ as shown in Fig. 1 from an 
explicit calculation. Since TlZ(a, 0) + a = A 2 — uo 2 , this 
minimum corresponds to a maximum value u — uo^ be- 
yond which no real solution for lo exists. 

It is interesting to note that these two values a* and 
are actually identical because their parent equations 
TdR/ da + 1 = and 1 — Tl = are so in view of the re- 
lation dR/da = —X. The significance of this coincidence 
is that the imaginary part of the self energy grows con- 
tinuously, so that there is no discontinuity in the density 
of states, and therefore when the band gap vanishes the 
resulting metal-insulator transition is also continuous. 

Solving for real oj from Eq. (33) for a given a > a* 
and (3 = 0, we thus have the solution to the Eqs. (33), 
(34). If to* is the solution corresponding to a = a*, then 
since a* = , we have lo* — , beyond which no real 
solution for lo exists, as discussed earlier. Therefore lo* 
is an uppcrbound for solutions with f3 = 0, for which 
the imaginary part of the self energy and the density 
of states vanish. This implies that lo* is the band-edge 
energy, and lo < lo* is the energy gap region, with the 
energy gap E g (T) = 2uj*(T). 
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FIG. 2. The normalized energy gap E g /U in two dimen- 
sions vs. the fluctuation strength in the range < Sn 2 < 1/2 
for different U values. The critical U value, above which the 
band gap never closes, is slightly above 10. 

A convenient way to obtain u*, and hence the band 
gap E g (T) is from the plot of TTZ(a, 0) + a vs. a in the 
range < a < A 2 . The variation of the band gap so 
obtained with the fluctuation magnitude Sn in the range 
< Sn 2 < 1/2 is shown for two and three dimensions 
in Figs. 2 and 3 respectively for different values of U. 
There exists a critical interaction strength U c , such that 
for U < U c with increasing Sn the energy gap vanishes 
at some critical value Sn c (U), whereas for U > U c the 
band gap remains finite and the system remains insulat- 
ing. In two and three dimensions we obtain the critical 
interaction strengths U c — 10.2 and 12.5 respectively, 
(U c /W = 1.27 and 1.04). 

With increasing fluctuation strength V the solution a* 
of the equation 1 = Tl(a, 0) increases, as also the right- 
hand side of Eq. (33), so that the corresponding to* de- 
creases until at some critical value T = T* we finally 
have lo* — 0, and the left-hand side of Eq. (33) cannot 
increase any further. This critical fluctuation strength 
at which lo* and hence the band gap vanish marks the 
metal-insulator transition. For r > T* a solution of Eq. 
(33) exists only if /? ^ 0, which leads to a finite density 
of states at lo = 0. It is important to note that (3 and 
the density of states increases continuously with T — T*, 
signifying a continuous metal-insulator transition. 

While the range a > a* corresponds to the energy-gap 
region with (3 = 0, the other side a < a* necessarily goes 
with (3 7^ 0, and therefore corresponds to the band region 
wherein the self energy is complex. It is also clear from 
Eq. (33) that for lo 2 > A 2 a solution exists only if a is 
negative, so that the band region corresponds to a < 0. 
To solve Eqs. (33), (34) for a, (3 for some fixed T, we first 
solve for j3 from Eq. (34) for a given a, and then solve 
for a real lo from Eq. (33) with this pair a, (3. 



FIG. 3. The normalized energy gap E g /U in three dimen- 
sions vs. the fluctuation strength in the range < Sn 2 < 1/2 
for different U values. The critical U value, above which the 
band gap never closes, is slightly above 12. 

By varying a between a sufficiently negative value and 
a*, the full range of solutions are thus obtained for lo > 
lo*. From these self-consistent solutions for a, (3 the com- 
plex self energy is then obtained from (A 2 — lo 2 ) + X(u>) = 
a + i/3. 



A. Density of states 

We now discuss the evaluation of the density of states 
using this complex self energy E(w). For this purpose we 
rewrite Eq. (28) for the Green's function as 



G(k,w) = 



LO - 



£ 2 (k,L 



£ 



£ 



£ 



£ 



(35) 



where £ 2 (k,io) = E 2 + E(w). Writing £(k,co) = 
|£(k,w)|e 4e ( k ><"\ from the solutions for \£(k,Lo)\ and 
0(k,w) 



\£(k,L0)\ 

6(k,Lo) 



[{£ 2 +£ r H} 2 + £ 2 ( W 
i 



il/4 



1 

- tan 
2 



E 2 + E r (w) 



(36) 



we obtain £(k, lo) = |£(k, Lo)\e 



i0(k,w) 



£ r (k, lo) 



i£i(k,Lo), in terms of which the density of states is fi- 
nally obtained from 



N{lo) 



-E 

2tt V 

k 



Si 



£i 



(Lo-£ r ) 2 +£ 2 (Lo + £ r ) 2 +£ 2 



(37) 



Here we have made use of the transformation prop- 
erty of ek which changes sign under the transformation 
k — > k + 7r while £(k, lo) does not, so that when summed 
over all k the £k term in Eq. (35) yields a vanishing 
contribution. 
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FIG. 4. The electronic density of states N(u>) for different 
values of the fluctuation strength Sn. 

The results for the density of states for a simple cu- 
bic lattice and U = W/2 = 6 are shown in Fig. 4 for 
several values of the fluctuation strength Sn. Due to the 
fluctuation-induced states in the gap, the band gap de- 
creases with Sn, and at a critical value Sn c w 0.4, the 
band gap vanishes and the system undergoes a continu- 
ous metal-insulator transition. 



VII. BROKEN-SYMMETRY STATE 

So far we have been concerned with the spin-symmetric 
state. However, the present approach provides for a de- 
scription of the broken-symmetry state as well, and in 
this section we briefly make contact with the known re- 
sults for the AF state with respect to the AF band ener- 
gies and amplitudes. 

We consider the NN hopping case, and proceed with 
the decoupling approach A. In the AF state the density 
(riii) nas va l ues ( n ~ m)/2 and (n + m)/2 on the two 
sublattice sites, where m is the sublattice magnetization. 
Translational symmetry therefore exists only within the 
sublattice basis. After employing the same decoupling 
procedure as before leading to Eq. (15), Fourier trans- 
formation within the sublattice basis leads to 



G(kw) 



1 



uJ 



mU 
2 



mU 
2 



(38) 



The quasiparticle band energies in the AF state are thus 
given by uj — ±i?k, same as for the symmetric case. To 
make contact with the HF results for the AF state, we 
make the following approximation consistent with the HF 
scheme 



[mi - 2 J ) « {mi g ) 



m 



(39) 



In this case A 2 = U 2 ((n tl - n/2) 2 ) = U 2 m 2 /4, so that 
2A = mU, which now involves the sublattice magnetiza- 
tion. With this approximation, Eq. (38) simplifies to 



G(kw) = 



1 



uj 2 -E 2 



u-A 



<~k 
-A 



(40) 



with i?k = 1/A 2 + e 2 ,. This is precisely of the same 
form as obtained in the HF description of the AF state. 
]l7t Consequently, the expressions for the AF-state elec- 
tron amplitudes are identical, resulting in the same self- 
consistency condition 1/U = J2w l/^kj which deter- 
mines the sublattice magnetization m. 



VIII. SPIN CORRELATIONS AND THE 
ELECTRONIC SPECTRUM 

AF correlations have already been considered while in- 
troducing the decoupling schemes A, B, and C in sec- 
tion IV. Here we examine the influence of short-range 
AF ordering on such aspects of the electronic spectrum 
as the charge gap, the density of states etc. In the 
strong correlation limit and for half-filling (n = 1), when 
charge fluctuations are absent, the expectation value 
(nil + n i+s.i — n) vanishes only for perfect AF order- 
ing when (ni+si) = (n^). However, for T > Tn, when 
the spin ordering is not perfect, the finite spin correlation 
length £ results in a twisting of neighbouring spins by 7r/£ 
on the average. If at site i + S the spin is pointing up, but 
at site i the spin instead of pointing down is twisted by a 
small angle 9, then (mi) — cos 2 (#/2) w 1 — 8 2 /4, so that 
(nil + n i+5l) is slightly reduced from 1. Whereas if the 
spin at site i is pointing down, but at the neighboring 
site i + S the spin instead of pointing up is twisted by 
the small angle 9, then (rijj_ +n i+ si) is slightly increased 
from 1. This results in a twist-induced fluctuation term 
arising from the expectation value (n^ + ni + $i — n) which 
fluctuates between —9 2 /4 and 6* 2 /4, so that 



\Sn^\ = \ (n ti +n t+ si - n)\ 



4£ 2 



(41) 



This twist-induced fluctuation term couples with the 
hopping term at first order in Eq. (14), so that the quasi- 
particle band energies are modified by an order-t term. 
This feature is similar to that in the spiral state of the 
doped Hubbard model, where the twisting of neighbour- 
ing spins activates the hopping term at first order and 
broadens the bands symmetrically. [[35-37| The pulled- 
down states are filled whereas the pushed-up states are 
empty due to doping, resulting in stabilization of the spi- 
ral state. 

This twist-induced fluctuation term will also con- 
tribute to the overall fluctuation strength discussed in 
section V, with the disorder strength T now including 
a term related to the inverse spin-correlation length. 
Therefore with increasing temperature the enhanced 
magnon excitation and the spin-fluctuation induced spin 
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disordering will also contribute to states in the gap and 
the band-gap reduction. Thus in general, both charge 
and spin fluctuations contribute towards an effective dis- 
order strength T which will determine the overall elec- 
tronic spectral properties. 



IX. DOPING AND TRANSFER OF SPECTRAL 
WEIGHT 



In the strong coupling limit, doping with holes or elec- 
trons results in transfer of spectral weight between the 
two Hubbard bands. Consider the atomic limit (t = 0) 
where this feature is most pronounced and easily under- 
stood. A site with a single electron of spin a contributes 
two states at energies e CT = (occupied) and = U (un- 
occupied) . When a spin-cr hole is introduced by removing 
the spin-cr electron, the spectrum changes to both states 
at energy e a = e— = 0, because electrons of either spin 
can now go into the empty site with zero energy. Extend- 
ing to a collection of sites, if x is the hole concentration, 
then the density of states (both spins) at energies and 
U are 1 + x and 1 — x respectively, so that the density of 
unoccupied states in the lower energy level is 2x. Finite 
hopping introduces a spread in the hole energy of order 
t, so that this picture of the transfer of spectral weight 
remains valid as long as the band separation (<~ U) is 
much bigger than the bandwidth (<~ t). 

In order to examine the effect of hole doping on 
the spectral weights, we consider Eq. (14) appropriate 
for the NN hopping case, and focus on the correlation 
function term -ill 2 (^\T{n ti (t) - n/2} 2 a^(t)a^(0)\^), 
which reduces to the single-particle Green's function in 
the half-filled case n = 1. For a total particle den- 
sity n = (riii) + {riii) = 1 — x < 1) corresponding 
to a hole concentration x, the term (n^ — n/2) 2 re- 
duces to (n/2) 2 + xnii, which now contains an opera- 
tor term for finite doping. However, the doping-induced 
term— xiU 2 (^Tn^ (^0^(^)0^(0)1*) can be written in 
terms of G(t) from Eq. (4) as 

-xiU 2 (*\T nil (t)a n (t)a^(0)\*) = 
xU[id t G tl {t) - 5(t)6ij + J2 tsGi+s^t)} , (42) 

8 

so that the correlation function term under consideration 
can be exactly obtained in terms of the single-particle 
Green's function G. Decoupling the equation for G as 
before in section IV A, we obtain after Fourier transfor- 
mation, 

[(w - nU/2) 2 - {(nU/2) 2 + e£} - xU(u> - e k )]G(kw) 
= lu- nU/2 + e k - xU . (43) 



G(kw) 



xU/2 + e k 



After simplification and the Hartree shift uj — U/2 
we obtain: 



where as before 



(U/2f 



xUeu 



(44) 



The quasi- 



particle band energies are therefore given by ±£ k = 
±y / £'^ — xt7e k , which in the strong coupling limit reduce 
to ±[£7/2 — xe k ], when only terms upto first order in t/U 
are retained. This indicates that doping leads to a qual- 
itative change in the nature of the band itself, with the 
dispersion changing from e|/(7 of order J to xe k , result- 
ing in an effective doping-induced hopping strength and 
bandwidth of order xt. This is again very different from 
the Hubbard I "rigid-band" result wherein there is only 
a slight doping-induced modification of the bandwidth, 
but no qualitative change in the nature of the band. 

We examine the loss of integrated spectral weight in 
the upper band due to doping. Integrating the spectral 
function A(u>) over the upper band, and summing over 
both spins, we obtain 

/e /-e i 

<1ujA{uj) = 2 J d^-Y^ ImG(kw) 



E 



g k - xU/2 + e k 



(45) 



Evaluating the momentum sum in the strong coupling 
limit by retaining terms only upto order t 2 /U 2 , we ob- 
tain, 



/■© 

! / dujA( 



t z 



x 2 ) 



(46) 



Upto order t the integrated spectral weight in the upper 
band is therefore simply reduced to 1 — a; by doping, as in 
the atomic limit. Similarly the weight in the lower band 
is enhanced to 1 + x. 



X. CONCLUSIONS 

In conclusion, the new decoupling scheme is appropri- 
ate for both the spin-symmetric (paramagnetic metallic 
and insulating) phases as well as the broken-symmetry 
AF phase. Although restricted to bipartite lattices, the 
decoupling approximations involve only weakly fluctuat- 
ing operators, resulting in several improvements over the 
Hubbard I scheme. Independent of magnetic ordering, 
the scheme yields, at the lowest order, the correct strong- 
coupling bandwidth of order J, a non-zero critical inter- 
action strength (above which the band gap opens) only 
if same-sublattice hopping (e.g., NNN hopping) is also 
present, and the correct quasiparticlc spectral weights 
(of 1/2 per spin) in each band for the half-filled model. 

Charge and spin fluctuations, which are manifested in 
double occupancy and vacancy of sites and spin twist- 
ing due to short-range AF ordering, were studied within 
a static, random approximation. The self-consistent, 
CPA-level evaluation of the disorder-averaged self energy 
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showed a fluctuation-induced band-broadening and con- 
sequent decrease in the band gap. The critical interaction 
strength, above which the band gap remains finite even 
when the fluctuations are included to their full extent, 
was obtained as U c /W — 1.27 and 1.04 in two and three 
dimensions respectively. 

The solution of the non-linear self-consistent equation 
for the disorder self energy shows that with increasing 
fluctuation strength the band gap shrinks to zero con- 
tinuously and the density of states N(0) between the 
bands also grows continuously. The monotonic depen- 
dence of fluctuation strength on temperature translates 
into a continuous metal-insulator transition. The fluc- 
tuation magnitude as well as its dynamics can be de- 
termined from the appropriate correlation function, and 
this is presently under investigation. The dynamics is 
relevant in view of the switching of the spectral function 
between the two Hubbard bands, and the possibility of a 
motionally narrowed central peak in the density of states 
when the Hubbard switching rate becomes comparable 
with the energy separation between the bands. Quan- 
tum Monte Carlo simulations do show a central peak in 
the vicinity of the metal-insulator transition. [ fi"H 

With doping a quantitatively correct transfer of spec- 
tral weight between the two Hubbard bands in the strong- 
coupling limit was obtained at the lowest-order level in 
the decoupling scheme. Furthermore, in contrast to the 
Hubbard I approximation, a qualitative change in the 
nature of the band itself was obtained, with an effective 
doping- induced hopping strength and bandwidth of order 
xt. 

Extension of the equation of motion approach to the 
transverse spin propagator to obtain the low-energy spin- 
fluctuation spectrum will provide a magnetic description 
of the spin-symmetric state, and enable one to approach 
the magnetic phase boundary from the paramagnetic 
side. This, as well as correlations and dynamics of fluctu- 
ations, and their impact on the electronic and magnetic 
spectrum, are presently under investigation. 
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